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Abstract 

We compute the order af (next-to-leading) corrections to the Wilson coeffi- 
cients of the AF = 2 effective Hamiltonian in the Minimal Super symmetric Stan- 
dard Model with completely arbitrary soft terms. These results are relevant for 
phenomenological studies of neutral meson mixings in the presence of large squark 
mass splittings, in particular in hierarchical scenarios of supersymmetry. These cor- 
rections achieve a considerable reduction of the uncertainty related to the matching 
scale ambiguity, and are therefore numerically relevant. We briefly analyze the ef- 
fect of certain mass splittings on the size of the NLO corrections, and compare it 
with the case of almost degenerate squarks in the mass insertion approximation. 

1 Introduction 

There is an almost absolute consensus among the particle physics community that un- 
known physics must exist at energies of around the TeV. Indeed, the minimal Higgs 
mechanism, so far completely untested, is quite unlikely to be the true mechanism for 
electroweak symmetry breaking; and even if it were, an additional explanation should be 
provided for the stabilization of the electroweak scale. Another indication of new physics 
around these scales is given by, for example, the dark matter problem. Low-energy su- 
persymmetry is a leading candidate to fill that gap since it provides natural and simple 
explanations to those conundrums. 

However, over the last decade, two other things have become clear. First, that the 
structure of flavor symmetry breaking of the Standard Model is a very peculiar one, and 
second, that, quite unexpectedly, flavor physics experiments confirm it to great accuracy. 
Indeed, any new physics model one can think of leads naturally to violent flavor breaking, 
and the experimental limits on flavor-changing neutral currents (FCNC's) are among the 
most stringent constraints that any model must satisfy. 

The Minimal Supersymmetric Standard Model (MSSM) is no exception: an enormous 
set of new flavor violating parameters arise from its supersymmetry-breaking sector, 
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which in the quark sector can be understood as a general misalignment between quark 
and squark mass matrices. More specifically, a simultaneous flavor rotation of quark 
and squark fields diagonalizing the Yukawa matrices leads to the so called super-CKM 
basis, in which no tree-level FCNC couplings are present. But in this basis the squark 
mass matrix M 2 is not generally diagonal, and an extra rotation of the squark fields 
is necessary to diagonalize this matrix, introducing strong interaction flavor violating 
couplings. Therefore, SUSY contributions to flavor violating observables mediated by 
strong interactions will easily compete with the SM contributions, which are driven by 
weak interactions, and well tested experimentally. Flavor constraints on these flavor 
violating parameters have been studied extensively, and the conclusion to be taken from 
those studies is, basically, that it is quite difficult to reconcile the relatively low SUSY 
masses required by naturalness with a generic flavor structure of the soft breaking terms. 

A solution to this problem is to assume that the structure of flavor violation of the 
new physics respects the hypothesis of Minimal Flavor Violation (MFV), by virtue of 
which the Yukawa couplings are the only source of flavor symmetry breaking. However, 
this is quite a pessimistic scenario for new physics searches in the quark sector, and it is 
useful for many purposes to go beyond it. Indeed, present indications of new physics in 
B ds and K processes [HOE] seem to require a departure from MFV. 

In the MSSM, a suitable approach for phenomenological studies that takes into con- 
sideration all these issues, is to assume that squark masses are very nearly degenerate. 
This case is parametrized by a squark mass matrix of the form (M 2 )ij = M 2 (l + 8)ij, 
where M s is an "average" squark mass, and 5 is a matrix with entries much smaller than 
one. These small parameters, or mass insertions, parametrize the departure from MFV, 
and constraints on their values derived from flavor observables provide valuable infor- 
mation for model building. They are usually treated as expansion parameters, which 
at leading order (linear and quadratic in S for AF = 1 and AF = 2 processes respec- 
tively), define what is coventionally know as the Mass Insertion Approximation (MIA) 
[I]. Strong bounds have been derived for these quantities from FCNC processes [51 El [7], [8] 
and vacuum stability requirements [9l|10] as well as from charged-current processes 

However, it is also of interest to consider other scenarios with non-degenerate squark 
masses. For example, it has been argued that a "hierarchical" setup in which the first 
two generations of squarks are much heavier than the rest of the SUSY spectrum (lying 
near the electroweak scale), can satisfy naturalness criteria [12]. In this framework, 
correlation patterns between AF = 1 and AF = 2 observables can be quite different 
from the ones in the degenerate case. Indeed, as shown in Ref. pj3], in these scenarios 
the bounds from B — > can be partially evaded, allowing for a large B s mixing phase. 
This would be of great relevance if the experimental indication for such a large phase 
[H [151 [2] is confirmed. 

In order to perform phenomenological studies of these scenarios, corresponding cal- 
culations of flavor violating processes have to be performed. Here we focus on AF = 2 
processes, that is, K — K, D — D and B d <s — B^ s mixing. These low energy observ- 
ables are more conveniently computed in the framework of an effective theory in which 
heavy modes have been integrated out. The most general effective hamiltonian relevant 
for AF = 2 processes is given in Eq. ([6]), where the Wilson coefficients Cj contain the 
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information from heavy modes. Then the observables are expressed as functions of the 
Wilson coefficients and the matrix elements of the operators, and arise from the AF = 2 
amplitude in the effective theory: 

A cS = J> (O t ) = £ C, (<% + ^n 3 + OK 2 )) (^)(°) , (1) 

i ij 

where we have written the matrix elements of the operators in terms of tree level matrix 
elements (O)( \ The matrix elements must be computed using some non-perturbative 
approach, for example in the lattice. The SUSY contributions are encoded inside the 
Wilson coefficients, which are evaluated by matching the full theory (MSSM) onto the 
effective theory at some matching scale \l. According to the renormalization group (RG) 
prescription for the resummation of large logarithms, the matching scale /i must be close 
to the SUSY scale, and the scale at which matrix elements are computed must be close to 
a relevant mass scale in the effective theory (for example mj in the case of B — B mixing) , 
and the Wilson coefficients at the matching scale must be used as initial conditions for 
the RG evolution that provides the Wilson coefficients at the low scale. This evolution 
is governed in particular by the anomalous dimensions of the operators. 

Leading order (LO) strong interaction matching conditions in the MSSM have been 
computed in Refs. [T6l [T7J [18], and arise from the squark-gluino box diagrams shown in 
Appendix lA.il The corresponding next-to-leading order (NLO) corrections arise from the 
two loop diagrams shown in Appendices IA.2I - IA.41 and have been computed in Ref. [TH] 
within the Mass Insertion Approximation. The anomalous dimension matrix for the 
complete set of operators (Eq. (JTj)) has been computed at NLO in QCD in Refs. [20| [21] . 

The purpose of this paper is to present the computation of the NLO matching con- 
ditions beyond the MIA, in the presence of arbitrary soft terms, and in particular for 
arbitrary squark mass splittings. The motivation for a NLO determination of the match- 
ing conditions is three- fold. First, LO matching conditions are both scale- and scheme- 
independent, so in order to get scheme independent results and NLO scale invariance it 
is necessary to go beyond the leading order. Second, the leading order corrections are 
proportional to a 2 s . Since at LO neither the renormalization scale or the scheme can be 
specified for the strong coupling, LO results show a particularly high uncertainty related 
to the scheme and scale ambiguities. Third, this uncertainty is particularly severe due to 
the large anomalous dimensions of the AF = 2 operators involved. These uncertainties 
are largely cured by the NLO corrections, from about 10-15% to a few percent, as shown 
explicitly in Ref. |19j . 

The AF = 2 amplitude in the MSSM up to NLO can be written as 

A MSSM = £ a* (*f + + <^> (0) , (2) 

% 

where and F^ are the LO and NLO contributions, and we have factored out a 
common a 2 s . The matching of the full theory onto the effective theory is performed by 
imposing that the effective and MSSM amplitudes are equal at and below a matching 
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scale (jl. An order by order identification of Eqs. ((Tj) and (jSJ) leads to the following formula 
for the Wilson coefficients: 

C^aiFf^^-^Ffr. + O^). (3) 

3 

Thus, the NLO matching calculation requires the computation of the matrix r and the 
functions Fj. The matrix r is obtained from the renormalization of the operators of 
the effective theory, and its computation is described in Section [2j This matrix has 
been computed before in several renormalization schemes (see Refs. [20 | [2TI [19] ). and we 
agree with their results. The functions F{ are obtained computing the one- and two-loop 
diagrams in the MSSM (see Appendix EJ) . The details of this computation are described 
in Section [3], and the functions obtained here are the main new results of this paper. 

Both calculations, the effective and the full theory amplitudes, are carried out in 
the NDR scheme (dimensional regularization with anticommuting 75), with modified 
minimal subtraction (MS) of ultraviolet divergencies. Also, we choose massless external 
quarks with zero external momenta. This choice of external states introduces infrared 
(IR) divergencies from diagrams in which a gluon connects two external legs. We regu- 
larize these divergencies with an unphysical gluon mass A. While both r and depend 
on the external states (and are therefore IR divergent), this dependence cancels in the 
Wilson coefficients, as expected. 

After describing the computation of the relevant amplitudes, we discuss briefly in 
Section H] some issues related to the renormalization of ultraviolet divergencies and the 
renormalization scale dependence. In Section [5] we mention some of the checks that can 
be done to ensure the correctness of the results. In Section [5] we show how reduce the 
exact results to the MIA and to the MIA with non-degenerate squarks, which allows 
to compare our results with those in Ref. [19]. Finally, some results are presented in 
Section 

Before getting down to brass tacks, we would like to set some notation and specify 
some definitions related with rotation and mass matrices beyond leading order. At tree 
level, the super-CKM basis is defined by doing a joint rotation in flavor space of quark 
and squark fields such as to diagonalize the tree level Yukawa matrices. The resulting 
squark mass matrix is not diagonal and defines the tree level mass insertions. In this basis 
there are no tree level FCNC's. This matrix is diagonalized by an additional rotation of 
the squark fields, 

a i,L Dl 3 ' a i,R ~ Dr J 

= T ul^3 > = T U*R^3 ( 4 ) 

where {q L , gjj) denote the squark fields in the (tree level) super-CKM basis, and q is the 
mass eigenbasis. The rotation matrices Tl and Tr are 3x6 matrices, and the indices 
U, D will be omitted hereon, which raises no confusion. In the mass eigenbasis, tree 
level FCNC's appear; for example a flavor changing qrqj-g vertex is generated with the 
following Feynman rule: 

-zg s V2T a (riP L -T%P R ) , (5) 
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where T a are the color matrices and Pl,r are the chiral projectors. 

At NLO, a subtlety arises because a squark-gluino loop generates a finite flavor- 
changing self energy for the quark fields (see Fig. [2]). At this point, one must specify what 
is meant by the super-CKM basis at NLO, since the definitions for the mass insertions 
depend on that choice. We believe that the most natural definition for the super-CKM 
basis is the one for which quark fields do not mix at one loop, and tree-level FCNC's 
are absent. In this case the rotation matrices, the mass insertions and the CKM matrix 
differ from the tree level ones. This criterium does not coincide with that in Ref. [11] . 
but it has no effect when comparing our results with Ref. [THj, since this is not an issue 
in the degenerate case. A comment in favor of the criterium adopted in Ref. pj] is that 
mass insertions are directly related with SUSY-breaking parameters in the lagrangian. 
A further discussion on this issue is provided in Section HI and the formulae necessary 
to switch from one criterium to the other is provided in Appendix O 

2 Effective Hamiltonian for AF = 2 processes at NLO 

The most general effective Hamiltonian for AF = 2 processes up to operators of dimen- 
sion six can be written as 

t=i i=i 

where Ci are the Wilson coefficients and Oi are the dimension six AF = 2 operators. In 
four dimensions there are eight independent operators of this type. Here we choose the 
following basis: 

01 = s al ^P L b a s^P L bp 

2 = s a P L b a spP L bp 

3 = s a P L b p S(iP L b a 

04 = S a P L b a SpP R bj3 

5 = s a P L b p s p P R b a (7) 

where Pl,r = (lT7s)/2 are the usual chiral projectors. The operators 0i,2,3 are obtained 
from 01,2,3 by exchanging L «-> R. To simplify the notation throughout the paper we 
focus on the case of B s — B s mixing; for the cases of K, D and mixing one should 
make obvious replacements of the quark fields. 

In order to define the effective Hamiltonian beyond leading order, one must specify a 
renormalization scheme. Here we choose to regularize ultraviolet divergencies in dimen- 
sional regularization, where d = 4 identities between operators do not hold. This means 
that one must complete the effective Hamiltonian with a set of evanescent operators that 
vanish in 4 dimensions but can give finite contributions beyond leading order if they are 
multiplied by a divergence. The choice of a set of evanescent operators is not unique, 
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Figure 1: One loop diagrams contributing to the matrix elements of the four-fermion operators 
in the effective theory. 



and different sets lead to different subtractions, so specifying this set is necessary to fix 
the renormalization scheme. Here we choose the set of evanescent operators given in 
Ref. [21J. By requiring that the matrix elements of evanescent operators vanish in four 
dimensions one can omit these operators altogether from the effective Hamiltonian once 
the renormalization has been performed. 

The calculation of the NLO effective Hamiltonian amounts basically to the NLO 
renormalization of the operators. The matrix elements of the bare operators and of the 
renormalized operators are related through the renormalization constants: 



(Oi) 



bare 



(8) 



where Z s ^ are the quark wave function renormalization factors and Zij is the renormal- 
ization matrix necessary to renormalize properly the operators in the effective theory. 
The effective theory amplitude can then be written as 



(9) 



'■j 



5Zlj and the second term contains the counterterms. 



where we have written Z'^ = 5{j 
In order to obtain the NLO renormalized operators (Oj) ren in terms of tree level matrix 
elements (Oi)^°\ one must compute the one loop gluonic corrections such as those shown 
in Fig. [TJ For UV divergencies we use naive dimensional regularization with modified 
minimal substraction (MS-NDR). Moreover, we choose zero external momenta and set 
the quark masses to zero. This simplifies the computation but introduces IR divergencies 
from the soft gluon region. We regularize this divergencies using a gluon mass A. The 
same IR divergencies should appear in the full theory, and cancel in the matching, 
providing a non-trivial check of the calculation. 

All the one loop diagrams contributing to the NLO renormalized operators are pro- 
portional to the same loop integral, and the amplitude can be written as 



A 



cS 



On , 

3 16tt V e 



21og(A//x) )A ij + 6Z' i 



(Oi) 



(o) 



(10) 



The computation then provides the coefficients Ay, which must be extracted up to 0(e) 



A 



ij 



A?- + eA\y The sum over j runs over both physical and evanescent operators, the 



later giving finite contributions to A, 
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The counterterms can now be fixed according to the MS scheme: SZ'^ = (a s /167re)A?-, 
which provides the leading order anomalous dimension matrix of the operators (JTj) in 
QCD: 



7 = ^ 



-i 



dZ 



a. (o) = a, 

47T 47T 



4C F 



2 . 



(11) 



dlogfi 

where we have used the well known QCD quark wave function renormalization factor 
Z q = 1 — (a s /47re) Cp. We obtain 

A o o o\ 

-28/3 4/3 

16/3 32/3 
-16 

-6 2 
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(0) 



(12) 



in agreement with, for example, refs. [26, 19J. The missing 3x3 block corresponding to 
the operators (5i,2,3 has been omitted: they do not mix with the other operators and 
their anomalous dimensions are the same as for (^1,2,3 ■ 

Finally, the NLO amplitude in the effective theory is given by 



(13) 
(14) 



For the NLO matrix r we obtain 



/ 


-4/3 








\ 




f-5 











\ 







-44/3 4/3 













1/3 


-1 













16/3 16/3 








log(A/ y u)+ 





-15/2 


-25/6 
















-64/3 
















19/3 


-3 


V 








-6 - 


-10/3 ) 




V 








-1/2 


-7/6 J 



The first term is the IR divergent piece that must cancel in the matching procedure. 
The second term is the NLO contribution, and it is scheme dependent: it is valid only 
in the MS-NDR scheme. However it can be easily translated to dimensional reduction 
(DRED) and Regularization-Independent (RI) schemes by using the formulae in Refs. 
[j~9| [20] . In fact, as pointed out in Ref. [20], the matrix r can be thought of as defining 
the renormalization scheme, and Eq. ([TBI as a definition the renormalized operators. 
The same scheme must be used in the full theory calculation and in the evaluation of the 
matrix elements in order to obtain scheme-independent results for physical amplitudes. 
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3 Details of the calculation in the MSSM 



The NLO AF = 2 amplitude in the MSSM is obtained by computing the one- and two- 
loop Feyman diagrams listed in Appendix [A] The one-loop contributions are the box 
diagrams shown in appendix IA.1( the computation of these graphs is standard and will 
not be reviewed here. Renormalization of the two-loop diagrams requires the inclusion 
of these boxes with vertex counterterms, as well as pentagons containing the self-energy 
counterterms. These one loop integrals must then be computed up to and including 
terms of 0(e), which provide the finite scale- and scheme-dependent contributions to the 
NLO amplitude. We will come back to this in section HI 

The computation of the two-loop diagrams is done in four steps: 

1. Partial fractioning of the denominators and tensor reduction to reduce the Feynman 
integrals to a set of scalar integrals with three scalar propagators. 

2. Decomposition of the scalar integrals down to a set of one- and two-loop mas- 
ter integrals. These master integrals are known functions of masses and contain 
divergencies up to 0{l/e 2 ). 

3. Cancellation of UV divergencies. 

4. Fierz rearrangement and spinor transpositions in order to express the spinor struc- 
tures in terms of tree level matrix elements of the operators. 



Let us describe in some detail each of these steps. 



1. Since we are choosing massless quarks with zero external momenta, there are 
only two independent momenta appearing in the Feynman integrals. One can then 
decompose the denominators such that each diagram can be expressed as a sum of terms 
of the following type 

f( M (- r' Mi ' Vi '-h \ <- r^ 2 ' V2 '-h \ f ^ D( ?i d°q2 g^gf 2 ■ ■ • q^q^ 2 ; ; ; 

nrns ){S/3 L a ^ o a)[ s y i a ^ o 5) j ^tt) d {qf - m^q 2 - ml)^(Aq 2 - m 2 3 )^ 

where the /'s are some functions that depend in general on all the masses appearing in 
the Feynman diagram (including the fictitious gluon mass), and Aq = q\—q 2 . Also, V and 
r represent Dirac and color structures. The spinors s and b might be u or v spinors and 
they might appear transposed and in different order, according to the chosen reference 
order adopted to keep track of the relative signs of interfering Feynman graphs |22j. 

The tensor integrals are momentum-independent (again, because external momenta 
are zero), so they can be expressed in terms of scalar integrals multiplied by metric 
tensors. For example, 

d D qi d D q 2 q»q\ _ 9^ f d D q x d D q 2 q x • q 2 
(2u) D (2u) D (■ ■ ■) D J (2ir) D (2tt) d (• ■ ■) 



S 



In this way we can express each diagram as a sum of terms of the form 

d D qi d D q 2 {ql) a {ql) b (qi-q2Y 



/(m's) (s/s r^'-; b a ) (s 7 rg*;;;. &*) ^ 



(27r) fl (27r) D (gf - m\) n ^(ql - ml) n '^(Aq 2 - mf)^ 

At this point, care must be taken when manipulating the Dirac structures V and T 
after contraction with the metric tensors. In particular, in D = 4 — 2e dimensions, 
structures of the type (7 M 7„ 7a • • ■) <£> (7 M 7 I/ 7 A • • •) cannot be reduced like in D = 4, and 
evanescent structures (of order 0(e)) must be introduced. According to the NDR pre- 
scription here adopted, however, one can freely anti-commute the 75 and use the usual 
anti-commutation relations for 7 matrices. 

2. The resulting scalar integrals can be further reduced to a set of one- and two- 
loop master integrals by the method of recurrence relations [23]. This reduction can be 
performed automatically using the Mathematica program TARCER [23]. In this way 
the scalar integrals can be expressed in terms of various one-loop tadpole integrals and 
a single two-loop master integral, 

z(m " m2 ' ma) = J bT^M^lW^l) • (15) 

The result for this master integral with arbitrary masses is given in ref. |25j . 

3. Once all the loop integrations have been performed, the resulting expression for 
the Feynman diagram consists of a sum of terms of the form 

A(mg,m) B(m~ g ,m) \ - 

f 2 + 1 + C ( m ~9> m ) J ( S P T aJ,... U ( S 7 T aJ,... b S) ( 16 ) 

where A, B and C are some functions of gluino and squark masses. Since in this case 
one-loop corrections are finite, no 1/e 2 divergencies can appear, so all 1/e 2 terms should 
(and do) cancel directly. We also get automatic cancellation of 1/e terms for the dia- 
grams shown in appendix IA.41 as it should be. The rest of the diagrams (those shown 
in appendix IA.3I) contain 1/e divergencies that cancel against one- loop diagrams with 
vertex and self-energy counterterms. We will discuss the details of the renormalization 
in section HI 

4. After all the divergencies have been removed, 4D Fierz identities can be used 
and transposition of spinors can be performed to put all spinor and Dirac structures 
in suitable form. These structures must then appear in the precise combinations that 
constitute the tree level matrix elements of the physical operators, as for example, 

2 «7m^<) itl.PLul) + 2 (u^P L vl) (v^PluZ) — > (OjW . 

This, however, does not occur for individual diagrams, but only for certain groups of 
diagrams (and of course for the amplitude as a whole). Therefore this step provides an 
interesting check of the calculation, in particular of the relative sign between the different 
diagrams. The NLO quantities F^ 1 ' in eq. (j2j) are then obtained by summing all the 
contributions. 
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4 Renormalization 



In order to reduce the ultraviolet divergencies arising from the two-loop graphs, appro- 
priate counterterms must be introduced. In particular, the squark-quark-gluino vertex, 
as well as the quark, squark and gluino wave functions and masses have to be renormal- 
ized. For reference, the relevant renormalization factors in the MS-NDR scheme up to 
0(a s ) are 

f> 1 rv 1 

An e 9 Air e 

Z §a = l-^l(2N c + C F ) , Z~ q = 1 + 0(a 2 s ) , Zq = l-^-2C F , (17) 
47r e An e 

which are defined in the usual way. Also, a non-diagonal squark-mass counterterm is 
necessary, because at one loop divergent flavor- changing squark propagators are gener- 
ated. There are two ways of dealing with this issue: 1) Renormalize the squark rotation 
matrices so that non-diagonal squark masses are zero at one loop. In this case all the 
finite pieces of the one loop corrections cancel with the renormalized non-diagonal mass 
insertions, which are not zero anymore but of order 0(az s ). 2) Renormalize the non- 
diagonal squark-mass parameters minimally, and include the finite pieces of the loop 
corrections to the flavor-changing squark propagator. In this case the squark rotation 
matrices are defined such that renormalized mass insertions are zero at the matching 
scale. 

We choose the second option, and renormalize the squark mass parameters according 

tO m 2 b are = m | + w ^ 



2C F 

Aire 



k,q,q' 



"^'(/"Matching) = . (18) 

A consequence of all this is that, even though we are working in the basis of diagonal 
squark masses, non-diagonal mass parameters do run with the scale and contribute to 
the renormalization group equation. 

There is an additional issue related with the fact that NDR is a regularization scheme 
that breaks supersymmetry. In particular, the coupling g s appearing in the quark-gluon- 
gluon vertex and the coupling g s that appears in the quark-squark-gluino vertex, receive 
different radiative corrections in the NDR scheme. In the effective theory only g s appears, 
so when doing the matching it is convenient to have the MSSM amplitude expressed solely 
in terms of this coupling. Being this is an 0(a s ) effect, one can set g s = g s in the two- 
loop amplitude, but it gives a finite contribution of 0(a 2 ) from the LO amplitude. At 
the end, this effect is corrected for by performing in the one loop amplitude the following 
replacement (see [2~7] [TU]). 

~ndr = ^ndrA + ^4\ (1q) 



10 




q' 



Figure 2: Flavor changing quark self energies mediated by a squark-gluino loop. 

Finally, one must address the wave-function renormalization of the external states. 
There are two types of contributions: gluon corrections and squark-gluino corrections, 
both contributing a factor of —Cf to Z q (see eq. (|17p ). In addition, squark-gluino 
corrections contribute finite pieces of two types: 

1. Flavor diagonal corrections: 

These give a finite contribution to the on-shell quark wave-function renormalization 



constant, 



log(m?/^ 2 ) - V f(m k /m~ g ) 



^kq* j^kq 
/ „ J ""gj >- L,R L.R 

k 



(20) 



with f(x) = (x 2 - Ax + 3 - 2x(x - 2) loga;)/2(a;- l) 2 . 

2. Flavor changing corrections: 

The presence of quark flavor changing squark-gluino loops (see Fig. [2]) can be 
handled in two different ways, in relation to two different definitions of the super- 
CKM basis: 

1) Tree-level definition of the super-CKM basis: This is the usual definition, in 
which the tree level Yukawa matrices are diagonal. At one loop, quarks of different 
flavors mix through the radiative corrections in Fig. [2j and one must include these 
corrections in the external legs. 

2) "On-shell" definition of the super-CKM basis: In this case the quark superfields 
are (finitely) renormalized with matrix-valued renormalization factors, that induce 
a perturbative rotation in flavor space. These counterterms are defined at each 
order in perturbation theory such as to render the quark self-energies diagonal at 
the given order, effectively canceling the flavor- changing self-energies of Fig. |2j 

As mentioned in the introduction, the mass insertions (and the squark rotation ma- 
trices) that arise in each scheme represent different quantities. Here we assume an 
on-shell definition of the super-CKM basis, and omit the flavor-changing external- 
leg corrections in our computations. For completeness, we provide in Appendix [C] 
the explicit relationship between squark rotation matrices and mass insertions in 
both definitions. The full results for the NLO Wilson coefficients consistent with 
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the usual tree-level definition of the super-CKM basis can be obtained by adding 
an 0(a s ) correction to the rotation matrices T l [ R in the LO results, as explained 
in detail in Appendix 

Once all the renormalization factors have been defined, one can write down explicitly 
the renormalization group equation for the Wilson coefficients, 







+ 



d a q 



d dm\ 
+ 



d log jj? d log jj? da s d log /i 2 dm 



8 dm% d 



d log jj? dm 



^2 



C7(/0 = 0. (21) 



The renormalization group functions are given by 



dmfj 



d log /i 2 



dm- 
d log /i 2 



m ; 



4tt " 9 



(2n f - 6AQ 



(22) 



k,q,q' 



d log /i 2 

Therefore, with the notation of Eq. ([2"lh the explicit form of the NLO RG-equation reads 







,(i) 



9 log /i 2 ' 



where appropriate sums over i,j, k, q, q' are understood. 



d 



(23) 



,(o) 



(0) c,(0) 



7,2 ' / h 2 7i ' ? 



(9m 



5 Checks of the calculation 

There are several checks that can be made of the NLO calculation: 

1. Cancellation of UV divergencies: All the UV divergencies arising from the two- loop 
graphs (specifically those in appendix IA.3j) must be cancelled by one loop graphs 
with the insertions of the counterterms specified by the renormalization factors in 
eqs. 071) and (jlgj). 

2. Projection onto tree-level matrix elements: The NLO amplitude must be ex- 
pressible in terms of tree-level matrix elements of the physical operators, that 
is, *A NLO = F^ 1 ' (Oi}W . As mentioned before, this does not happen for individual 
two-loop graphs, and provides a check of interference between diagrams. 

3. Cancellation of IR divergences: Both the amplitude in the effective theory and the 
amplitude in the MSSM depend on log A, where A is the gluon mass introduced 
to regularize the IR divergencies. This dependence must cancel completely in the 
matching when combining both amplitudes; that 

is, FP - if >r« must be IR- 
finite. This is a non trivial check involving three completely independent pieces of 
the calculation: the matrix r from the renormalization of the effective operators, 
from the one-loop graphs, and f!> from the two-loop graphs. 
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4. NLO renormalization scale independence: The RGE (eq. (1231) ) must be fullfilled. 
There is a close relationship between cancellation of UV divergencies and fulfillment 
of the RGE, but it is nevertheless a convenient check. 

We have verified that our results fulfill these requirements. Moreover, a final check 
consists in reducing the results to the degenerate MIA and comparing them with the 
results obtained in ref . [19] . In the next section we address the issue of how to reduce 
the exact results to the MIA and the NDMIA. 

6 Reduction to MIA and NDMIA 

6.1 Reduction to the Mass Insertion Approximation 

Any Wilson coefficient computed in this paper has the structure of a sum of terms of 
the following type, 

f^(rh 2 , m), ml) ■ T% Pj • rg* rg • r£* Tf 

/ 6 , 2 (mf , rh 2 , ml) ■ T? rg ■ rg • T k * (24) 

f (rh 2 rh 2 rri 2 rri 2 \ ■ V iq * V ib ■ V js * V jq ■ V ks * V kq> ■ V lq '* V lh 
f a n(m 2 m 2 m 2 m 2 \ ■ Y iq * Y ih ■ T js * T jb ■ F fcs * V kq ' ■ V lq '* V lq 

where A, B, C, . . . are either L or R, and a sum is understood running over the indices 
i,j,k,l = dL,SL,bL,d R ,SR,bn, and q,q' = d,s,b, (for / 6j2 also q = u,c,t and k = 
u L ,c L ,t L ,u R ,c R ,t R ). 

These terms can be translated into functions of the entries of the squark mass matrix 
in the super-CKM basis (M 2 ) using the following relations, 

E r T T b = ; £ mf r f T% = [(M 2 )"]^ . (25) 

i i 

In the MIA, the diagonal elements in M 2 are assumed to be equal to an "average" 
squark mass M 2 , and the off-diagonal elements (called mass insertions and denoted by 
A^f), are assumed to be much smaller than M 2 . In this way, a power expansion on the 
dimensionless mass insertions = A^?/M 2 1 can be made, keeping in this case 
only the leading terms 5 2 . The mass eigenvalues are then M 2 (l + 0(5)), so the functions 
in Eq. ( |24l can be Taylor-expanded around the average squark mass: 

f(ml rh 2 , ...) = f(M 2 , M 2 , . . .) + f^-\M 2 , M 2 , . . .) • (rh 2 — M 2 ) + - - • 

+ fM-\M 2 , M 2 S , ...)■ (rh 2 - M 2 s )(m 2 - M 2 S ) + ■ • • (26) 



13 



Using the relations (1251) it is easy to see that, for example, 

h^mlrhl ml) • V% ■ V°* r? • T% = (27) 

tiT 1] {MlMlM* s ) A™ A™ + f^\MlMlM*)A L *A% , 

and similar for the other terms. 

We have checked that our results, when reduced to the MIA, reproduce exactly the 
results in Ref. [T9] . 



6.2 Reduction to the Non-degenerate MIA 

From the exact results obtained in this paper, one of the most straightforward gener- 
alizations of the MIA results that can be obtained is the case in which the diagonal 
elements of the squark mass matrix are non-degenerate. In this case, squark masses can 
be widely different while keeping mass insertions small. We call this the N on- degenerate 
Mass Insertion Approximation (NDMIA), in which diagonal masses are non-degenerate 
but mass insertions are kept only up to second order. 

In order to reduce the terms in Eq. (1241) to the NDMIA, we first expand the functions 
around m = and apply the relations ( 1251) . 

. , . , . „ f (ni,n 2 , ■■■)( n\ 

f(ml ml ...)■ vr r7 ■ rr rj? • • • = E „ ,„ , li^rd m 2 Y%i ■ ■ • , 

ni,n 2) ... 

and then we split the mass matrix in two pieces: M 2 = X + A, where X is diagonal and 
A <C X. Then we can resum the series in X, and keep only terms quadratic in A. For 
example, 

f(ml ml ml) • rf r» • rf r£ • r*r rf = f^(x h , x~ bR , x~ SR , x h ) A™A? b L , 

with the function F given by 

F s 2 J ( x s L > x i R ; x s R , x i L ) = 

/ ( X s L > X s L ; ) ~ / > ; X sr ) ~ / ( X S L ; ^ > ) + / > > ^ ) 

(A Si - X- bR )(Xs R - XgJ 

From the NDMIA results one can also recover the Wilson coefficients in the MIA by tak- 
ing the limit in which all X^s are equal. In the non-degenerate case, the dimensionless 
mass insertions can be defined normalizing by any of the different squark mass parame- 
ters. It is customary to normalize by some "average" squark mass, although it makes no 
real difference as long as in the degenerate limit the definition coincides with the one in 
the MIA. We shall understand that they are normalized by the smaller diagonal entry 
of the squark mass matrix. 
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Figure 3: Relative importance of the full NLO result with respect to the LO, for C% and 
C3, as a function of Xh = rh^/m-, where m 12 is a common mass for first and second 
generation squarks. The dashed lines correspond to the degenerate (MIA) scenario. We 
have chosen nig = ji = 350 GeV, and 400 GeV for third generation squark masses. 



7 Results 

The full results for the LO and NLO Wilson coefficients in the MS-NDR scheme are 
presented in Appendix [B] In this section we present some results for the NLO Wilson 
coefficients in a scenario with a hierarchy of masses, making a comparison with the MIA. 
In line with the rest of the paper, we keep focusing on the B s system for illustration. 

The considered scenario is obtained from the mass insertion approximation with 
non-degenerate diagonal entries, taking a common mass for the first two generation 
squarks (rn.12), different from a common mass for third generation squarks (777,3). This 
corresponds to the "hierarchical" scenario of Ref . [33] , where 771 3 is assumed to be near 
the electroweak scale, and m\2 is allowed to be heavy, up to several TeV. We therefore 
denote Xh = fh\ 2 /rrv~ and x\ = rhl/m^, for heavy and light respectively. Note that 7fi 12 
and m 3 are related to the true masses by corrections of 0(5). In this case the mass 
insertions are normalized to 7ri 3 , which corresponds to the MIA definition when Xh — > x\. 

The plots in Fig. [3] illustrate the relative importance of the NLO corrections, as a 
function of the mass splitting between the light and heavy squarks. The NLO correction 
is typically a ~ 10% effect in the degenerate case, but its importance increases with 
the mass splitting. For heavy squarks of about a TeV, the NLO contribution to C3 
can be up to a ~ 25% correction. However it should be mentioned that this is true in 
the NDR scheme and could vary in other schemes. Also, this depends on the matching 
scale fi, but a high sensitivity to fi would be related to large scale ambiguity at LO 
that is efficiently reduced at NLO. This means that the NLO contributions are either 
numerically important, or they achieve a considerable reduction of theoretical errors, 
and both situations are not easily disentangled. 

In order to analyze more closely the role of the mass splittings in the NLO corrections, 
we consider Ci, C 4 and C 5 as a function xi, for different splittings between xi and x^. This 
is shown in Fig. HJ where the dashed lines correspond to xi = Xh (that is, the degenerate 
case), and -departing smoothly from that limit- the solid lines show increasing values 
of Xh/xi. In these plots we take 5l,r = Sr,l = and m g = \i = 350 GeV. We 
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1.5 2.0 2.5 3.0 1.5 2.0 2.5 3.0 1.5 2.0 2.5 3.0 

Xi Xi Xi 

Figure 4: NLO Wilson coefficients Cf LO , Cf LO and C^ LO as a function of x\ = rn|/m~, 
where m 3 is a common mass for third generation squarks. The plots are in units of 
{o.1/^)8 2 LL and {0.1/11)8^8^ for C\ and C 4j5 respectively. The different lines correspond 
to different values of x^ = m^/mj, and range from the MIA (dashed) to 

x h = l.bxi, 3^, 4xi, 5xi, 6x t . 



see that increasing the heavy scale tends to reduce systematically the size of the NLO 
contribution, being largest in the degenerate case. 

As mentioned before, this is scheme dependent (although the conclusion might be 
more general). In any case, the true impact of the NLO corrections can only be estab- 
lished by analyzing their effect on observables. A full phenomenological analysis of these 
corrections and their impact on the bounds on the mass insertions (beyond the mass 
insertion approximation) is worthwhile, and will be presented elsewhere. 



8 Conclusions 

In this paper we have presented the computation of the NLO strong interaction correc- 
tions to the Wilson coefficients relevant for AF = 2 processes in the MSSM, beyond the 
mass insertion approximation. The full results for the Wilson coefficients in the NDR 
scheme are given in Appendix [B] 

These results are relevant for two reasons. First, NLO corrections are necessary to 
cancel renormalization scheme and scale dependence from the renormalized operators. 
This has the effect of a considerable reduction in the theoretical error. Second, in order 
to study scenarios with significant mass splittings one must depart from the degenerate 
mass insertion approximation. We have shown some illustrative examples of the effect 
of mass splittings in the NLO Wilson coefficients as compared to the degenerate case. 

A full phenomenological study of neutral meson mixing incorporating these new cor- 
rections will be presented in the future. Corresponding calculations for AF = 1 processes 
are underway, and when available will allow to perform a complete analysis of correlations 
between decay and mixing observables of neutral mesons at NLO in a s . 
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A List of diagrams 

In this appendix we present the set of one- and two-loop feynman diagrams that con- 
tribute to the MSSM amplitude. We focus on the case of B s mixing; for other neutral 
meson systems the external quarks must be changed. Dashed lines and i,j,k,l denote 
squarks, solid arrowed lines and q denote quarks and solid lines without arrows are 
(majorana) gluinos. We show all the possible topologies with the different insertions 
of external quark flavors, but an additional multiplicity is present interchanging initial 
and final external quarks (whenever they correspond to different contractions). For ex- 
ample, at leading order there are two diagrams of each of the two topologies shown in 
Appendix IA.ll with appropriate signs given by a reference order of external legs. 

A.l LO diagrams 



A. 2 IR divergent diagrams 
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A. 3 UV divergent diagrams 
A. 3.1 Vertex corrections 
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A. 3. 3 Squark self-energies 
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A. 4 Finite NLO diagrams 
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B Wilson coefficients 



In this appendix we present the complete expressions for the Wilson coefficients at NLO 
in the MS-NDR scheme. We omit the results for Ci,2,3j which are obtained from Cx,2,3 
by exchanging L «-> R in the squark rotation matrices. The coefficients are split into LO 
and NLO: 

a(/i) = C?V) + CfV), (28) 

and depend on the matching scale \i through a s , m- g and rhi, and explicitly through 
log(m?//i 2 ). The functions depend on the squark masses through the variables Xi defined 
as %i = mf/mj. 

The LO Wilson coefficients are given by 



Cf(A*) 



at 



12 2 m? 



(s G lj - 22 H ij ^j r£*r£ r j L s *r j * 



at 



68 G ij vrrv,' r^Trf 

12 2 m? 



a 



1 O -pis*-pib -pjs*-pjb 

122m? G L * L i> L * Ll 



at 



12 2 m? 



168 G tJ + 24 iP J ' J rj?Tg rfTj + 44 ^ T l L s T^ r£*r£ 



az 



12 2 m? 



8 G ij - 40 H ij ) T%*T% T j L s *T^ + 60 H ij rjTg r£T£ 6 



(29) 



where a sum is understood over 2, j = d L , sl, b L , cIr, sr, 6r. The functions G, are: 



G l 



2xi log Xi — xf + 1 

\Xi 1) (Xi X 'jj 



+ (Xi <-> x,) ; H" 



2x 2 i log Xi + 3x 2 — 4xi + 1 



+ (Xi <-> Xj) 



The NLO coefficient functions are the main result of this paper. They are scheme 
dependent; we present the results in the MS-NDR scheme, but can be translated easily 
to other schemes using formulae analogous to that presented in Ref |19j . The function 
Li 2 (x) denotes the dilogarithm defined in the usual way: 



Li 2 (x) = - 



dt 



log(l - 1) 



t 



(30) 



Appropriate sums over indices are understood, in particular q, q' = d,s,b and Q = 
u, d, s, c, b, t, and i,j, k run over left and right-handed squarks of u or d-type depending 
on the quark they appear with in the rotation matrix: for example a term containing 



20 



contains a sum over k = UL,CL,tL,UR,CR,tR. All dependence on w-type quarks and 
squarks come only from diagrams with gluino self-energies. The NLO Wilson coefficients 
read 



12 3 7rm? 



^kb 
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_i_ ijk pjs*pj& pis*pj<7 -pfc(?*p/cfe _|_ Z,*.?^ pis*pi& p?s*pj'& /pfcQ*pZsQ . p/cQ*pZsQ\ 
yi,2 L L L R R L ' 1 L L L L \ L L R R ) 



-m-i 
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+$} (rf rj? rj*rg p^Tf + (l <-> #)) + (rf rg rg*rg r^Tf + (l^r)) 

i iijfc pis*pi& pjs*pjfr /p/cQ*p&Q , p£;Q*p/cQ\ , i ij'/e p^s*pi& pjs*pj& /p^Q*-pA;Q , pfcQ*pfcQ\ 
^ 4,1 Z, L R R \ L L ' R R / ' 4,2 R R L L \ L L ' R R J 

+h% (Pf p* rj*ri 6 (r£ Q *p£ Q + vf*vf) + {l^r)) 
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The coefficients g ljk , h ljk , m ljk , n ljkl and p ljkl are the following, 



tf* = 15 (Af + Af ) - 25 Af + UAf + (j <-> fc) 

#f fe = -8(Af + ) + 81 (Af + Af ) + 34Af - 25 Af + (j <-> fc) 

#f = 34(Af + Af - 2Af - 2Af + Af ' + Af ) - 133(Af + Af ) 

#f = -6(Af + Af - 2Af ' - 2Af + Af + Af ) + 15(Af + Af ) 

sf = -ll(Af + Af + 7Af + 7Af ) - 1 33 (Af + Af ) 

g f 2 = -84( Af + Af ) - 21Af - 6Af + (j <-> fc) 

$J = 36(Af + Af ) - 6Af - 21Af + (j <-> fc) 

sf = -ll(Af + Af + 7Af + 7Af ) - 133(Af + Af ) 

yg'f = -15( Af + Af + 7Af + 7Af ) + 15( Af + Af ) 

^f = -4(Af + Af - 18Af - 18Af ) + 71Af - 62Af + (j <-> fc) 

<?g 3 fe = 12 (Af + Af ) - 62Af + 71 Af + (j <-> fc) 

^f = -l5(Af + Af + 7Af + 7Af ) + 15(Af + Af ) 



/if = -85(5f + B[ lk ) + ll(5f + 5f ) ; /if = -68(5f + B{ lk ) ; 

hf = 12(5f + B{ ik ) ; /if = -84(5f + B{ ik ) ; hf 2 = -12(5f + tff ) ; 

/if = 77(fif + B[ lk ) - U(B l 2 jk + B 3 2 lk ) ■ /if = -4(5f + B{ lk ) ; 

hl k = -UA(B[ jk + B{ ik ) + 20(£f + fif ) ; /if = 105(fif + B{ ik ) - 15(fif + B 3 2 ik ) ; 



m f = -2if f j + Cf ) + llF k {C l 2 j + Cf) ; mf = -17if + Cf) ; 

mf = 3if + Cf ) ; mf 2 = -llF k (C l 2 j + Of 5 = -15F fc (Cf + Cf ) ; 

mf = -21F fe ((f j + Cf) - 6if C^' + Cf) ; mf = -F k (C\ 3 + Cf) + lOif Cf + Cf) ; 



raf = -2{D\ jkt + D{ lkl ) + ll(L>f ' + Df L ) ; nf = -17(D\ jkt + D{ lkl ) ; 

nf l = 3(L>f + L>f ' ) ; nf = -ll(L>f ' + D j 2 lkl ) ; nf 1 = -15(D% kl + D j 2 ikl ) ; 

nf = -21(£>f + D{ ikl ) - 6(L>f ' + £>f ') ; ngf = "Pf + + 10pf ' + D j 2 ikl ) ; 
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p ijki = 4E u E ik . p m = _ 2E ij E ki . pijfci = 342£ lfe £ j7 ; 
p^' fc = -8E il E jk + l7lE ik E jl - E ij E kl ; pgf = 6£ ij £ H : 
p« fc = -8E il E jk - £ ifc £ j7 + 3^ J '£? W ; pg*' = -2E ik E jl ; 



with the functions A?'*, Cf , Dp', and F k given by 



2(xj - l) 2 (xj - l) 2 (x k - l)(xi - xj) 



(x 3 -l) 2 (x t -x k ) 2 U 2 (l-^) 



-8(xi - l) 2 (xj - l) 2 Li 2 (l - Xi) - 8(xj - l) 2 (x k - l)xi\ogXi 

+ (x k - l)(xi - x^XiXj - x,ix k - XjX k - Xi - Xj + 2rr fe )(4Li 2 (l - x k ) - log 2 x k ) 

+6(xj - l) 2 (xi - x k ) 2 \og 2 x k - A(xi - l)(xj - l)(x k - l)(xi - Xj) 

-4(x, - l)(xj - l)(xi - Xj)x k \ogx k - 8(xj - l) 2 x 2 logXilogx fc j 



A^ k 



- -f^ + ^-im-im-ix*-*,) t - a{x ' - 1)2(11 " 1)4 log2 x 

+(xi - l) 2 (xi - Xj)(x k - l)(2xj - x k - l)(4Li 2 (l - x k ) - log 2 x k ) 

-4(xj - l)(xi - Xj)(xl - l)x k \ogx k - 8(xj - l) 2 (2xfclogx fe - 5x k + 5)x- logx, 

-4(xj - l){xi - 6)(xj - l)(x k - l)(xi - Xj) 



1 1 9 

Aijk _ _ Aijk , _Aijk i f 

3 2 1 + 2 2 (xj - l) 2 (xj - l)(x fc - l)(xi - Xj) L 



(Xj - l)(x fe - l)Xj log 2 Xj 



+2(xj - l)(x fe logXfc - 2x fe + 2)xjlogXi + (xi - l)(x; - x-,)(x fe logx fe - 2x fc + 2) 



(xj - l) 2 (x i - l)(x fe - l)(x j - x fe ) 



8(x fc - l)(x; - x j ) 2 Li 2 (l 



+4(xi - l) 2 ^- - x fe )Li 2 (l 



X7. 



8( Xj - l)(x fe - l)(2xi - Xj - l)Li 2 (l 



+A(x k - l)(xi - Xj) 2 log 2 Xj - 8(x fc - l)xf logXj logXj 



+8(xi - l)(x k - l)xj logXj 
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B{ jk = 



r>ijk 
n 2 



jjijkl 



^2 



E ij 



(Xj — 1 — X; logXj)Xj logX 



2 4 (xi- l) 2 {xj - l)(xj - x k ) 
_%l A ijk + Mxi - l) 2 Li2(l - Xj) - (xj -x fc ) 2 log^ Xj + A(x k - 1) (a;, log a;, -Xj + 1) 



) 2 log 2 . 
2{xi - l) 2 {x k - 1) 



{Xi Xj) 



A(xi - x k ) 2 U 2 {l - g) + Xi(-4(2 + (-3 + Xi)xi)x k - Xi \og Xi 



-4 \ogXi(x k + Xi \ogx k ) + 21ogx fc (-2(2 + (-3 + x,i)Xi)x k + (x t - 2x k 

+ (3 + (-3 + Xi)xi)xl) logx k ) + A(xi - 2x k + (3 + (-3 + Xi)xi)x 2 k )U 2 (l - x k )) 



ijk , QxiX k (x k log x fc - 2 logx k + 2x fc Li 2 (l - x fc ) - 2) 



-( Xi -6)B?* + 
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2x; log Xi — x 2 + 1 



64x;(log(m?//i 2 ) + log^ - 1) 



C* j = -16 



2a^ log Xi — 3xf + 4xi — 1 



{Xi 1) (Xi Xj) 

2xi(x k - l) 3 logx; 



(Xi - l) 2 (x k - l) 3 (x; - Xj)(x k - Xi) 

+ - l)((x k - l)(-x k (x k - 3) + x 2 (x k + 1) - Xi(3 + x 2 )) - 2(xi - l) 2 x fc logx fc ) 



32^(log(m?//x 2 )+log^-l) 



2x\[x k - l) 3 logXi 



(Xi - l) 2 (x k - lf(Xi - Xj)(x k - Xi) 

+(xi - l)((x k - l)(xi + x 2 (l - 3x k ) + 3xiX 2 k - x k (x k + 1)) + 2(xi - \) 2 x\ \ogx k ) 



2xi log Xi 



2xj \ogXj 



{Xi \){x,i Xj) \Xj \){Xi Xj) 



?k _ x\ - 4x fc + 3 - 2x k (x k - 2) \ogx k 
4(x fe -l) 2 



Finally, the functions f % * appearing in the part with four squark rotation matrices 
are given by: 



1 



(Xi - l) 3 (Xj - l) 3 (Xi - Xjf 



I 8(xi - Xj)(xj - 1) 3 ( — 115 - 72xi + 399x 2 - 2\2x\ 



+ (-32 - 220^ + 74x 2 + 88x 3 ) log x { ) log(m|//i 2 ) 



+ 
+ 



A(xi - Xj)xj(2395 + Xj(-3587 + 3579x; - (x t - 1)(1189 + llx^Xj + 3(x; - l) 2 x 2 )) 
(xj - l)(-3xf + xf(53 + 3Axj) + x 4 (700 - 1584^ + 773x 2 ) 
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+xl (102 - U38xj + 2107^ ~ H37a; 3 ) + a; 2 (-72 + 392^- + 422a; 2 + 231a; 3 - 34a^ 4 ) 
+x 2 (120 + 40xj + 50a: 2 - 3a; 3 ) + 2^(256 - 1082^- + 86a; 2 - 16a; 3 + 3a; 4 ))] log 2 ^ 



-2 



2(xi - l)(xj - Xj)(xj - l)(-983 + Xj(1232 + ^(419 + ^(-50 + 3^))) + 1890^ 



-^(2263 + ^(838 + 31xi))xj + (-1114 + 5^(274 + 73a;;))a; 2 ) Li 2 (l - x t ) 



+ 



2(xi - l)(xi - Xj)(xj - l)(-3a£ + a; 3 (50 + 31^) + ^(120 + Xj(40 + (50 - 3xj)xj)) 



+a; 2 (40 - BOxj + 94a; 2 ) + Xi(120 + a^-410 + Xj(-B0 + 31^)))) Li 2 (l 



+ - A(xi - Xj)(xj - l) 2 (-64(x j - 1) + ^(1563 + ^(-311 + ^(-119 + Sx t )) 

-1575xj + Xi(Ml + 98xi)xj + 3(x { - l)x 2 )) \ogXi 

+ (xi - l)(xj - l)(-3^ 5 + ^(50 + MX,) + x 3 (40 + a^-98 + 151^)) 
+x 2 (120 + Xj(40 + (50 - Zxj)xj)) + 2xiXj(32 + Xj(-81 + ^(-49 + 17 Xj))) 
+x 2 (120 + xj(-162 + Xj(-328 + 151^)))) log x< log xj } 

( 1,;,^" '',^,,, ..., r;) , { [ 8 0* " " X ) 3 (- 515 + 889 ^ - 437a; * + 63x * 

+2(-136 - 5xi + 39a; 2 ) logXi)] log(m?//i 2 ) 

+ A(xi - x j )x j (7A19 + ^(-12118 + 14097^ + (x { - l)(-4699 + 1979^)^)) 

+ 2( Xj - l)(-26l4 + x 4 (547 + 668^) + 2^(1088 - 3764a:,,- + 1470a; 2 + 261a; 3 ) 
+2x 3 (909 - 3013a;., + 979a; 2 ) + a; 2 (146 + 2862;, - 261a; 2 ) 
-2x 2 (610 - 809^ - 2156a; 2 + 1320a; 3 ))] log 2 ^ 



-2 



A(xi - l)(xi - Xj)(xj - 1)(— 918 + 261a; 3 + (1429 - 682a; j )a; j 



-llx?(26 + 37xj) + 2^(1325 + 38a; j (-47 + 28xj))) Li 2 (l - x t ) 



+ 



A(xi - l)(xi - Xj)(xj - l)(261x 3 - lla; 2 (26 + 37^-) 



+Xi(-U6 + (1156 - 407a; i )x j ) + 2^-146 + a^-286 + 261^-))) Li 2(l - ^) 
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+ 



+x;(-3823 + Xi(1129 + 518x; - 2165XJ-) + 4341^)) logx 



+ 



2(xi - l)(xj - l)(261xf + x 3 (-286 + 372^) + 2x^(139 + 6x^(25 + Zlxj)) 
-2x 2 (73 + 10x^-15 + Uxj)) + x 2 (-146 + a;, (-286 + 261^))) logx; logx^j 

(x , - mxj ~- mxi - Xl) > { i s{x > " - 1)3(121 - 187x - + 47l? + lte - 



+ (48 + 62xi - 74x 2 ) logx 4 ) log(m?//i 2 ) 



+ 
+ 



60(x; - Xj)xj(-47 + Xj(62 - A5xi + (x* - 1)(15 + 17^)^))] 
- 2(xj - l)(9x, 5 + x\ (-95 + 68xj) + x 2 (86 - 86^- + 9x 2 ) 



+2^ 3 (143 - 191xj + 57x 2 ) + x i x i (192 - 556^- + 346x 2 - 9x 3 ) 
-2x 2 (22 + YlZxj - 412x 2 + 208x 3 ))l log 2 x; 



A(xi - l)(xi - Xj)(xj - 1)(162 + 9x 3 + x 2 (-86 + 77 Xj ) + x^-247 + 94x^- 



+x;(-239 + 478 Xj - 248x 2 )) Li 2 (l - x*) 



+ 



4(xj - l)(xi - Xj)(xj - l)(xi + Xj)(86 + 9x- + x^-86 + 9xj) 



+x,(-86 + 68xj)) Li 2 (l - 5t) 



+ 
+ 



4(xj(557 + Xi(-395 + 222x* - 49xj) - 335 Xj ) - 96(x j - l))(x; - x j )(x i - l) 2 logx 
- 2(xi - l)(xj - l)(9xf + x 3 (-86 + 308x^0 - 2x 2 (-43 + 74x^- + 200x 2 ) 



+x 2 (86 + Xj(-86 + 9xj)) + 2x^(31 + 2xj(-37 + 77xj))) \og Xi logx^j 

{Xi - mx ~- mXt - Xj) 2 { [ ~ 48 ^ ~ x ^ ~ 1 ) 3 ( 167 - 333 ^ + 237x " ~ 71x " 
+ (112 - 91xi + 50x 2 + x, 3 ) logXi)] \og(ml/fi 2 ) 



+ 



48(xi - x j )x j (1354 + x^-2225 + 2616x; + 4(x* - 1)(— 217 + 97x i )x j + 3(x* - l) 2 x 2 )) 
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+ 



6(xj - l)(6xf + x\ (9 - 30^) + 4(121 - 254^ + 184a; 2 ) + 3x 2 (2 + 5x 2 + 2x 3 ) 



-2^(719 - 1324xj + 530x 2 + 132a; 3 ) - 3^(896 - 1750^ + 908x 2 + 45x 3 + 6xj) 
+2x 2 (339 + 2Q2xj - 1376x 2 + 898a; 3 + Ihxf))] log 2 ^ 



12(xi - - a; j )(a; j - l)(-620 + 6x- + 3a; 3 (5 - 8^) + (1219 - 572^)^ 



+x 2 (-50 + (103 - 26a; i )a; i ) + ^(613 + 2x j (-613 + 281^))) Li 2 (l - Xi 



+ 



36(xi - l)(xi - Xj){xj - l)(2x- + x 3 (5 - Sxj) + 2xj + x 3 A5 + 2xj) 



WiXjil + 8xj) + Xi(2 + Xj(-12 + Xj - 8a; 2 ))) Li 2 (l - 



+ 



- 24(xi - Xj)(xj - l) 2 (-224(a; j - 1) + 3^(1305 + ^(-610 + X;(35 + 6^)) 



-1248^- + 4^(127 + Xi ) Xj + 6fa - l)x 2 )) 



log Xi 



+ 



- 6(xi - l)(xj - l)(6a;f + 15a£(l - 2 Xj ) + 4a; 3 a; i (-31 + 16^) 



+3a; 2 (2 + x]{h + 2 Xj )) - 2x^(112 + £,,(-126 + a; i (62 + IhXj))) 
+2x 2 (3 + Xj-(126 + xji-77 + 32xj))))] log a^ log | 

4(xi - - 1) 3 (597 - 15333;; + 1539a; 2 - 603x 3 



-1 



(Xi - l) 3 (Xj - l) 3 fa - Xj)' 



+2^(352 - 509X; + 289a; 2 ) logo;*) log(m|//x 2 ) 



+ 
+ 



32(xi - Xj)xj(A87 + Xji-754 + 801 Xi + fa - l)(-267 + 47x i )x J -)) 



2fa - l)(27a; 3 (7 + 9xj) - 27xf (-9 + 25xj) - 27x^(5 + Al Xj + 34x 2 ) 



+x\ (-880 + 1409^ + 1631x 2 ) + a; 3 (65 + imixj - 5047a; 2 - 1035x 3 ) 
+53; 2 3; i (-143 + 367^- + 5053; 2 + 135rcf))l log 2 Xi 



Afa - l)fa - Xj)(xj - 1)(709 + 27s; 3 (-9 + 25xj) + ^(-1175 + 3Axj) 



+x 2 (871 - 3^(1445 + 722^)) + ^(-1067 + 3^(1405 + 958^))) Li 2 (l - Xi ) 



+ 



lOSfa - I) fa - Xj){xj - l){-{x 2 A7 + 9xj)) + x 3 (-9 + 25 Xj ) 
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+x 2 (-7 + (25 - 66 Xj )xj) + XiXj(-2 + 25^(1 + Xj ))) Li 2 (l - ^) 



+ 
+ 



8xi(xi - Xj)(xj - 1) 2 (1061 - 453 Xj + ^(-455 - 761 Xj + ^(274 + 334^))) log^ 
2(xi - - l)(27a£(9 - 25a;.,) + 27a£(7 + 9xj) - ^^(487 + 27^(34 + 25xj)) 



+xl (189 + xj(-918 + 377xj)) + x-Xj(-487 + 3^(2542 + 377^))) log a;; log a^j 

(^TH^W^) 1 { t 16(:ri " " 1)3( " 77 " 81Xt + 321x * " 163 ^ 



+ (-16 - 179^ + 58s; 2 + 65a;, 3 ) log a*) log^/^ 2 ) 



+ 



+ 



I6(xi - Xj)xj(1258 + Xj(-2015 + 2292a;; + 32( Xi - l)(-23 + 8^^)^ + 21 (a* - l) 2 ^))] 
- 2(xj - l)(42xf - x 4 5 (257 + 178^) + (-689 + 2414^- - 792x 2 ) 



+x 2 (-54 - 160^- - 215a; 2 + 42a; 3 ) + 2x 3 (119 + 292 Xj - 1182a; 2 + 660a; 3 ) 
+x i a; i (-128 + 634*;, + 780a; 2 + 37a; 3 - 42a;)) 

+2x 2 (21 - 410^ + 352a; 2 - 658a; 3 + 89a;))) log 2 Xi 

4( Xi - \){xi - X j)(xj - l)(-756 + ^(787 + a;,(178 + ^(-215 + 42a;,))) 

+ 15893^- - 3^(1382 + 3^(143 + 136 Xi )) Xj + 2(-610 + 3^(551 + 77a;i))x 2 ) Li 2 (l - x t ) 

+ A( Xi - l)( Xi - Xj)(xj - l)(-42x- + a; 3 (215 + 136xj) + a; 2 (160 + ^(-533 + Itexj)) 

+Xj(54 + Xj (l60 + (215 - 42xj)xj)) 

+ Xi (54 + Xj (-1Q4 + Xj (-533 + 136^))))] Li 2 (l - 



+ 



- 8(xi - Xj)(xj - l) 2 (-32(3; j - 1) + x;(1551 + 7x^-82 + Xi (-13 + 6xi)) 



-I728xj + 4(253 - 53xi)xiXj + 42( Xi - l)x 2 )) 



logs;; 



+ 



2(xi - l)(xj - l)(-42xf + x\ (215 + 178^) + 4a; 3 (40 + a^-183 + 32xj)) 



+a; 2 (54 + 3^(160 + (215 - 42a; j )x j )) + 2x 2 (27 + a^(-18 + 3^(59 + 64^-))) 
+2a;ja;j(16 + Xj(— 18 + Xj(— 366 + 89xj))))j logs;, loga;., j 
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+2xi(-32 + 7xi + 13x 2 ) logo;,) log(m?//x 2 ) 



+ 
+ 



96(xi - Xj)xj{-lh7 + Xj(214 - 171x* + (ar< - 1)(57 + 43x i )x j )) 
- 30(x i - l) 2 (9x 4 5 - 9x) + 9x,x 2 (3 + 2 Xj ) + 5x*(-16 + 7x i ) 



+xf(19 + 134xj - 63x 2 ) - x 2 x,,(65 + 16^- + 9x 2 )) 



log 



-2 



60(xi - l)(xi - x i )(x j - 1) 2 (— 71 + Xi(97 + £,(19 + 9x, - 46^; 



—70xj) + 62xj) Li 2 (l — Xj) 

+ [ - 540(x, - l) 2 (x, - Xj )\ Xj - l) 2 ] Li 2 (l - |) 

+ 24x;(x; - x i )(a: i - 1) 2 (515 - 435xj + x;(-137 - 23xj + X;(22 + 58^))) 
+ [ - 30(x; - l)(xj - l)(9x t 4 (x 3 - - 1) - 9(xj - l)x 3 + x 3 (3 + Xj-)(3 + 5xj) 



log Xj 



+x*x 2 (5 + 9xj(2 + xj)) + x 2 xj(5 + Xj(-7A + hx 5 ))) 



log Xi log X, 



C Flavor changing quark self-energies: Relationship 
between the tree-level and on-shell definitions of 
the super-CKM basis 

In this appendix we elaborate on the relationship between the two different definitions 
of the super-CKM basis addressed in Section HI 

The flavor changing self-energies of Fig. [2] can be written a£] 

S^(p) = ^ q (p 2 )PL + T$Hj?)Pr + j> l^, L q (p 2 )P L + X™(j> 2 )P R \. (36) 

When inserted as external legs, the quark propagator will provide a chiral enhancement 
with the chirality-flipping part of the self-energy, and we will keep only this contribution. 
Moreover we will expand in the external momentum, which is justified by the fact that 
the SUSY masses are much higher than the external quark masses. Then, the self energies 

1 These corrections can be found in Ref. [11]. A discussion of how these can be absorbed into wave- 
function counterterms is given in Ref. [28j . 
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we are interested in are given by 



= f^E^H) T n> k L q n (37) 

k 

with g(x) = (1 — x + xlogx)/(l — x). 

For the case of B s mixing, the corrected couplings of the external quarks are, 



ig s ^/2T^ a [WlP L -W^P R )u a h 



where, in the approximation m& ^> m s ^> m^, and keeping only the chirally enhanced 
terms, 

Ug = F id E LR /mb + F is E LR /mb (3g) 

nr = rf E^/m, - rf E^/m 6 (39) 

and correspondingly for L <-> i?. 

At this point, we denote by T^ ) the squark rotation matrices associated with the 
tree-level definition of the super-CKM basis, and T the rotation matrices in the on-shell 
case. At tree-level both coincide; at the one loop level, the relationship between both 
definitions is given by 

it = if* + it Tfflm + rt 

rf = ls * + rf z%/m s - rf e£?m (40) 

and correspondingly for L <-> i?. The full results for the NLO Wilson coefficients in 
terms of r^ ^ -corresponding to the tree-level definition of the super-CKM basis- can be 
obtained by making the substitutions of Eq. (l4Tj|) in the LO expressions given in Eq. (|29|) . 

In the degenerate MIA case, these corrections cancel. Consider for example the 
coefficient C\ in Eq. (1291) . Making the substitutions of Eq. ( 1401) in the rotation matrices 
with index i, and making the reduction to the MIA (see Section [671]) gives, at order af, 
two terms: 

/(*„ • • -)nr It = -^/(MfH, ■ ■ ■) g'(M>/ml) 8™ (41) 
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£ f( Xi , ■ • -)rr K = 2 -^f(M*/ml ■ ■ •) gXM*/ml) JfL 5™ (42) 

and both cancel due to the relative signs in Eqs. (j4T)l) . The same happens with the other 
two terms arising from the T J ^, and in the other Wilson coefficients. Therefore, the two 
bases coincide in the mass insertion approximation with degenerate squarks. 

Finally, the relationship between the mass insertions defined in each basis can be 
deduced from Eqs. ( |4"U1) and ( |25i) . As we have just seen, both are equivalent in the 
degenerate case. In the non-denerate case, one such relation is (at leading order in the 
mass insertion expansion), 

\LL a (0)ll 2a s m~ g g{X- 8L /m])-g{X- h Jm]) 

Asb ~ A ° b ~ 37^ X- SL - X- bR ^ ~ X '* ] Asb ' (43) 

with the notation of Section 16.21 Analogue relations for other mass insertions can be 
found accordingly. From this last relation one can also see that the effect disappears in 
the degenerate case. 

The effect of the squark-gluino corrections to the external legs (or the difference be- 
tween the two definitions of the super-CKM basis), turns out to be numerically important 
because of the chiral enhancement. While this article was under revision, a phenomeno- 
logical study of these corrections has been performed in Ref. [29], with the conclusion 
that the bounds on some mass insertions (corresponding to the tree-level definition of 
the super-CKM basis) are modified considerably. These corrections have obviously no 
effect on the phenomenological bounds derived for the mass insertions associated to the 
on-shell definition of the super-CKM basis, and in this scheme the NLO corrections are 
well defined in the limit of vanishing quark masses. However, it is not a claim of the 
present paper that the on-shell scheme is in any way preferred to the tree-level one. In 
fact, the tree level scheme might be more convenient if one wishes to relate phenomeno- 
logical bounds on these type of low energy processes to specific mechanisms of SUSY 
breaking. This issue has been already discussed in Ref. [29]. In those cases one can add 
the flavor- changing self-energies as given in that paper to our NLO results, or do the 
substitutions of Eqs. (jjOj) or (j4"3"]l directly. 
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